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ABSTRACT 

To understand the formation and evolution of solar-type stars and planets in the solar neighborhood, we need to obtain their stellar 
parameters with high precision. We present a catalog of precise stellar parameters for low-activity FGK single stars in a volume- 
limited sample followed by the HARPS spectrograph in the quest to identify extra-solar planets. The spectroscopic analysis was 
completed assuming LTE with a grid of Kurucz atmosphere models and using the ARES code to perform an automatic measurement 
of the line equivalent widths. The results are compared with different independent methods and also with other values found in the 
literature for common stars. Both comparisons are consistent and illustrate the homogeneity of the parameters derived by our team. 
The derived metallicities of this sample reveal a somewhat different distribution for the present planet hosts, but still indicates the 
already known higher frequency of planets observed for the more metal-rich stars. We combine the results derived in this sample with 
the one from the CORALIE survey to present the largest homogeneous spectroscopic study of the metallicity-giant-planet relation 
using a total of 1830 stars. 

Key words. Stars: fundamental parameters planetary systems Stars: abundances Stars: statistics 



1. Introduction 

More than half of a thousand planets have been detected since 
the discovery of t he fir st exoplanet orbiting a solar-type star 
dMavor & OuelozLll995h . The method that has been adopted the 
most widely in achieving these detections is the radial veloc- 
ity technique. There are several dedicated observing programs 
that are, almost, continuously monitor the sky. The data obtained 
by these programs consists of stellar spectra collected by differ- 
ent high resolution spectrographs located at many different ob- 
servatories around the world. Today the HARPS spectrograph 
(Mavor et al., 2003) is at the top of the "food-chain" when talk- 
ing about planet hunting. This is not only due to its high spectral 
resolution but more importantly its long-term stability in the de- 
termination of the radial velocity of stars. 

These new discoveries of planets are allowing us to con- 
strain the theory of formation and evolution of planetary 
systems. Soon after the first detections of planets, it was 
found that the stars hosting the newly discovered giant plan- 
ets were systematically more metallic than a normal sample 
of stars. This well-established correlation has been con firmed 
by several author s ([G onzalez, 1997; Gonzalez et al., 12001 



Santos et all l200 lL 12004: -Fischer & Valentk ,2005: IJdrv et al- 
2006t lUdrv & SantoT l2007n . and has been shown to imply 



that core accretion is the main mechanism for the formation 
of giant planets (I da & L in, 2004; Mordasini et al., 2009) and 
not disk-instabilitv dBossi [20021) . This observational correlation 
might also be consistent with more recently developed theories 



such as the form ation of planets by the tidal downsizing of giant 
planet embryos dNavakshinl |20 1 Ol) . 

In this paper, we present a catalog of spectroscopic stellar pa- 
rameters for FGK single stars of low activity that belong to the 
volume-limited sample observed with HARPS to search planets. 
In Sect. 2, we describe briefly the HARPS GTO sample ana- 
lyzed in this work. Section 3 describes the procedure used to de- 
rive precise spectroscopic stellar parameters as well as estimate 
the stellar masses. In Sect. 4, we compare our derived parame- 
ters with other independent methods described by our team in 
previous works. In Sect. 5, we compare our spectroscopic de- 
rived parameters with others that we found in the literature to 
check consistency between different methods, while In Sect. 6, 
we present the solar-twins that are present in the analyzed sam- 
ple. In Sect. 7, we combine this HARPS sample with the initial 
CORALIE sample. We compare both to check for differences in 
the [Fe/H] distributions. These two samples combined form the 
largest volume-limited sample of single stars of low activity with 
homogeneous and precise measurements of [Fe/H] that we use 
to study the metallicity correlation in planet host stars. Finally, 
in Sect. 8 we summarize the work presented here. 



2. The sample & observations 

The sample of stars presented here is part of a HARPS GTO 
program that aims to detect and obtain accurate orbital elements 
of Jupiter-mass planets in a well-defined volume of the solar 
neighborhood. Although giant planets were the main targets of 
this program, there is also the possibility of detecting Neptune- 
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Fig. 1. Distribution of the signal-to-noise ratio of the final spectra 
compilation. 



mass planets. This sample was compiled with stars in the so- 
lar neighborhood out to 57.5 pc ( ILo Curto et al.LlTOl O') from the 
Sun to complete a planet-se arch surve y previously started with 
the CORALIE spectrograph dUdrvet a l.. 2000). The selection of 
the stars within this limited volume (containing more than 1400 
stars) was caiTied out to ensure that the stars are suitable to ob- 
tain precise radial velocities. They have a wide range of spectral 
types extending from the cooler MO to the hotter F2, although 
the majority are solar- type stars. All the stars were selected to 
have low levels of activity and low rotation rates. The sample 
was also clean of known binaries and variable stars. The total 
number of stars within this program reaches ~ 850. 

This sample indeed represents an extension of the one stud- 
ied by the CORALIE program, and collectively they form 
the largest volume-limited sample with homogeneous measure- 
ments of [Fe/H]. In the past few years there have been some 
discoveries of planets orbiting the stars in this sa mple, includ- 



ing al so a few that were posteriorly unconfirmed ( Tinne v et al 



201C; 



20031; iHi brard et al., 2010; Santos et al., 2010; Naef et al 
Moutouetal., 2009; Lo Curto et al., 2010; Santos et al 
Mordasini et al., 2011; Mouto ueta l., 2011). These will allow 
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us to increase the number of stars with planets in a subsample 
to ensure a proper statistical comparison of the metallicity dis- 
tributions. 

The spectra for the sample were collected by HARPS during 
several runs, thus there was the need to combine the spectra of 
each star into one to achieve a higher signal-to-noise (hereafter 
S/N). The individual spectra of each star were reduced using the 
HARPS pipeline and combined using IRAF Q after correcting 
for their radial velocity. The final spectra have a resolution of R 
~ 110000 and values of S/N that vary from as low as ~ 20 to 
as high as ~ 820, depending on the amount and quality of the 
original spectra. 

Figure [1] shows the S/N distribution for the stars for which 
we were able to derive precise parameters. Although there is a 
significant amount of spectra have a low S/N, the vast majority 
(~ 75%) have a S/N greater than 100 and ~ 25% a very good 
S/N higher than 200. 

Our method for deriving spectroscopic stellar parameters is 
based on a differential analysis with respect to the Sun (using an 



adopted value of log^(Fe) - 7.47^3, therefore there is some lim- 
itation on the spectral type of the stars for which we can derive 
good and precise stellar parameters. Our spectroscopic method 
is effective for F, G, and K stars with temperatures ranging from 
~ 4500 K to ~ 6400 K. For this reason and since the sample is 
composed of a wider variety of spectral types, we made a pre- 
selection of the stars for which it would be appropriate to use 
our method to derive precise spectroscopic parameters. The se- 
lection was made by co nsidering the index co lor B-V taken from 
the Hipparcos catalog dvan LeeuwenL l2007h . We only selected 
the stars with B-V < 1.2 to avoid the coolest stars in the sample. 
This B-V cutoff reduced the number of stars to 65 1 . Amoung 
these 651 stars there are a significant number (69) for which we 
were unable to derive good and precise parameters. The main 
reason for this is that the spectra of most of these stars have 
a very low S/N. There were also some cases of stars with high 
rotation that affected the spectra in a way that prevented an accu- 
rate derivation of spectroscopic parameters using our automatic 
methods (e.g. problems with strong blended lines). Among these 
69 stars, there were also a few M stars that have passed the B-V 
cut-off and for which we cannot derive parameters. Altogether, 
the final number of stars with precise and homegeneously deter- 
mined parameters is 582. 

3. Stellar parameters 

The spectroscopic stellar parameters and metallicities were de- 
rived following the same procedure used in pre vious works 
(ISantos et alll2004ISousa et alll2006ll2008ll201 Ih . The method 
is based on the equivalent widths (hereafter EWs) of Fe I and 
Fe II weak lines, by imposing excitation and ionization equi- 
librium assuming LTE. We used the 2002 version of the code 
MOOG (iSnedenl [1973) and a grid of Kur ucz Adas 9 plane- 
parallel model atmospheres IKutucJ, I1993|) . In this procedure 
[Fe/H] is used as a proxy for the metallicity. The EWs were au- 
tomatically measured with the ARESQ code (Automatic Routine 
for line Equivalent widths in stellar Spectra - Sousa et al., 2007j) 
that reproduces with success the common "manual" and interac- 
tive deter mination of EW measurements. The procedure used in 
ISousa et al. (2008) was closely followed and the same input pa- 
rameters for ARES were used in this work. There is a significant 
number of stars in this sample with S/N lower than 100 and even 
a few spectra with S/N lower t han 50, to which w e applied the 
same procedure as discussed in'Sous a et alJ (1201 ih . 

There are 14 stars in common between this vol ume-Umited 
sampl e and the metal-poor sample analyzed in ISousa et aP 
(l20lTh . When comparing the results of these 14 stars, we ob- 
tained mean differences of -0.5 + 0.6 (cr = 8.6) K for effective 
temperature; 0.01 + 0.00 (cr - 0.02) dex for surface gravity; and 
0.00 + 0.00 (cr = 0.004) dex for [Fe/H]. This small mean dif- 
ferences demonstrates the high precision of our spectroscopic 
automatic analysis when using the same instrumental configura- 
tions. 

The eiTors were determined in the same way as in previous 
works. We recall th e discussion about p recision versus accuracy 
errors presented in ISousa et aO (l201 In . We present here preci- 
sion errors, although in the electronic table we also add the "ac- 
curacy" error 



' IRAF is distributed by National Optical Astronomy Observatories, 
operated by the Association of Universities for Research in Astronomy, 
Inc., under contract with the National Science Foundation, U.S.A. 



^ This value taken from iGonzalez & Laws! ( l200(]h was first used in 
Santos et al. ( 2000) and is adopted here to guarantee the homogeneity 
of our metallicities 

' The ARES code can be downloaded at 
|http://www.astro.up.pt/~sousasag/areS| 
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Fig. 2. In the top panel, we present the distribution of the sam- 
ple stars in the H-R diagram. We also plot some evolutionary 
ti-acks computed with CESAM for a 1.0, 1.1, and 1.2 Mq star 
of solar metallicity. The filled circles represent the planet hosts 
in the sample. In the bottom panel, we present the metallicity 
distribution of this sample. 



3.1. Masses and luminosities 

St ellar masses were e stimated following the same procedure as 
in ISousa et al ] (1201 Ih . where we applied the stellar evolution- 
ary models from the Padova group using the web interface deal- 
ing with stellar isochrones and their deri vatives H- We used th e 
Hipparcos parallaxes and V magnitud es (Ivan LeeuwenL l2007h . 
a bolometric correction from .Floweij (Il996h . and the effective 
temperature derived from the spectroscopic analysis. The errors 
presented for the masses were also acquired from the web inter- 
face. 

The luminosity was computed by considering the estimated 
Hipparcos parallaxes, V magnitude, and the bolometric correc- 
tion. Its error is derived based on the parallax errors, which are 
the main source of uncertainty in the calculation of luminosity. 
The typical eiTor in the luminosity is ~ 0.05, which was obtained 
assuming the mean parallax for the stars (~21 mas) and a typical 
error in the estimated parallaxes of 1 . 1 mas. 

Figure |2]presents some characteristics of the sample. The top 
plot shows the distribution of the sample on the Hertzsprung- 
Russell diagram, where we represent evolutionary tracks for 1 .0, 



1.1, and 1.2 M p stars of solar metallicity, compu ted using the 
CESAM code (iMoreli 1 1 997t iMarques et al.L 12008') R Moreover, 
we present the typical error boxes on this specific diagram, 
where the error in the luminosity was described just before and 
the error in the temperature is derived from our spectroscopic 
method. This plot shows that the sample is composed mainly 
of main-sequence solar-type stars but also contains some sub- 
giants. Two stars can be distinguished well below the main se- 
quence. These stars are identified as HD 84627 and HIP 99606, 
being the second the farthest from the main sequence. We be- 
lieve that the reason for these strange positions in the HR dia- 
gram comes directly from the parallax values that have indeed 
atypically large error values, ~ 40 + 7 mas and ~ 25 + 8 mas, 
respectively, for HD 84627 and HIP 99606. The Hipparcos vari- 
ability flag of HIP99606 is defined to be 3, meaning that variabil- 
ity is higher than 0.6 mag (the only one in the sample with flag 
set = 3). We indicate that HD84627 was found to have duplicity- 
induced variability in the Hipparcos catalog. 

In the bottom panel of Figure |2] we present the metallicity 
distribution, which has a mean value of about -0.10 dex. As 
expected, this value is compatible with both comparison sam- 
ples presented be fore in previous works (Santos et al.L I2004t 
ISousa et all l2008l) an d the original sample of CORALIE, as can 
be seen in the work of ISantos et al.l (l2004h and later in this work 
in the top panel of Fig. |7] 

In Table [T] we present a subsample of the derived spectro- 
scopic parameters. The full table will be available in electronic 
format. 



4. Internal comparisons 

To check the consistency of our determined parameters, we 
derived the effective temperature using two other independent 
methods developed in our previous works. 

We first used the latest calibration for the effec tive temper- 
ature as a function of B-V and [Fe/H] derived in 'Sousa et ajj 
(201 1). We used the color B-V from the Hipparcos catalogue and 
the spectroscopically derived [Fe/H]. The comparison between 
our spectroscopically determined effective temperature and the 
temperature derived using this calibration can be seen in the up- 
per plot of Fig. [3] The temperatures are quite consistent between 
them having a mean difference (Teff_Cal - Teff_Spec) of only 
-19 K with a standard deviation of 76 K. There is a slight off- 
set for the hotter stars, but this offset can easily be explained by 
the calibration that we use only being valid for effective tem- 
peratures ranging from 45QQK to 6400/r. This offset starts to 
be seen for stars with temperatures higher than 6100 K. When 
considering the stars with temperatures lower than 6100 K, we 
obtain a more relieble result with a mean difference of -5 K and 
a standard deviation of 68 K. In this comparison there is a clear 
outlier, which is identified as HIP99606 an outlier also in the HR 
diagram. There might be a problem with the color B-V obtained 
from the Hipparcos catalogue for this star. 

The second independent method that we used to check the 
consistency of our determined spectroscopic parameters was the 
one based on line ratios. We used ARES to derive the EW of the 
lines used to derive the line ratios. The estimation of the effective 
temperature was carried out using the freely available code "Line 
Ratio Calibration Code" accessible from the ARES webpag^ 
The method behind this code and the respective determination 



http://stev.oapd.inaf.it/cgi-bin/param 



' http://www.astro.up.pt/corot/models/ 
* The "Line Ratio Calibration Code" 
|http://www.astro.up.pt/~s ousasag/areS| 



can be downloaded at 
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Fig. 3. Top panel: Comparison between the spectroscopically de- 
rived effective temperature and the temperature derived with a 
calibration based on B- V and [Fe/H] . Bottom panel: Comparison 
between the spectroscopically derived effective temperature and 
the temperature derived using line ratio calibrations obtained 
with the "Line Ratio Calibration Code". 



of the cahbrations for each line ratio is described in detail in 
ISousa et aL There were some stars for which we could 

not derive the temperature using this method, for reasons similar 
to the one described before. There is also a range of temperatures 
that defines the validity of each line-ratio calibration. The typical 
range of temperatures is [4500, 6400] K. Nevertheless, we were 
able to derive the effective temperature for 564 out of 582 stars 
with an impressive mean difference of only 6 K and a standard 
deviation of 76 K. This not only shows that our measured stellar 
parameters are consistent and on the same temperature scale as 
our previous measurements but again proves that the line-ratio 
code is very efficient and precise because it derived consistent 
temperatures for a sample of FGK stars that were not used for 
the calibration of the fine ratios. 
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Fig. 4. Comparison between the spectroscopically derived effec- 
tive temperat ure (a and c) and [Fe/H] (b and d) and those pre - 
sented in No rdstrom et al I (I2004.) and lCasagrande et all (l201lh . 
respectively. 



5. Other studies 

We compared our results with those of other studies that have 
stars in common. This is useful to access the consistency be- 
tween the various methods that use either spectroscopy or 
photometry to d erive stellar parameters. We used the spectro - 
scopic results of lEdvardsson et all (ll993h . lBensbv et all (l2003h . 



and 'Valenti & Fischeil (|2005[). and the photometric results of 
Nordstrom et al. (2004) and lCasagrande et alJ (1201 ih . 

Figure |4] compares our measurements of the effective 
temperature (panel (a)) and [Fe/H] (panel (b)) with those 
of [N ordstrom et al. (2004), which used the calibrations of 
lAlonso et al.i(il996 ) to derive the effective temperatures. In these 
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shows a very low mean difference of only 13 + 7 (cr = 100) 
K and at higher temperatures a clear increase in the dispersion, 
which is expected as the stars are very different from the Sun 
an d the errors increase acc ordingly. The metallicities presented 
in ICasagrande et al.l (1201 ll) are compared in panel (d) from the 
same figure. The comparison is quite consistent with a mean dif- 
ference of only -0.01 + 0.01 (cr = 0.09) dex. 

Figure |5] shows a comparison between our different spectro- 
scopically derived parameters, effective temperature (top panel), 
surface gravity (middle panel), and [Fe/H] (lower panel) and 
other sp ectroscop i c val ues derived by other w orks, namely 
Edvardss on et al.l d 19931) (6 star s in common), Bensbv et al.l 
(2003) (3 stars in common), and IValenti & Fischeii (l2005h (28 
stars in common). 

For effective temperatures, we observe a mean difference of 
-41 + 11 (cr = 67) K meaning that our temperatures are on av- 
erage cooler. There seems to be a small offset for the hotter stars 
(with Teff > 6100 K) derived bv I Valenti & Fisched (l2005h . The 
reason for this offse t is unc lear. It is also present in a similar 
comparison made in ISousa e t al. (2008) but is not visible when 
comparing with the results of other methods considered in the 
same work. 

When comparing the surface gravity, we obtain a cloud of 
points around the identity line that is rather typical of surface 
gravities determined with different methods. The comparison 
finds a mean d ifference of -0.10 + 0.03 (cr = 0.18) dex. The val- 
ues derived bv Edvards son et aP (11993 ). which used the Balmer 
discontinuity index, seems to be sys tematic lower than our de- 
termination as already pointed out in lSousa et alj (l2008h . 

Finally, the comparison of [Fe/H] that can be seen in the 
lower plot of Fig.|5]shows a very good consistency with a mean 
difference of 0.00 + 0.01 (cr = 0.04) dex. This shows that the dif- 
ferent authors have very consistent metallicities measurements 
despite the differences in surface gravity and temperature, for 
the hotter stars. 

All these comparisons confirm the consistency of the param- 
eters derived by our team. 



Fig. 5. Comparison between our determined spectroscopic stel- 
lar parameters, effective temperature (top panel), surface grav- 
ity (middle panel), and [Fe/H] (bottom panel) and the mea- 
surements of other authors. The crosses, triangles, and dia- 
monds re present the compar i son p oi nts with the values deter- 
mined by 'Valenti & Fischer| (l2005h . iBensbv et all (l2003h . and 
lEdvardsson et al. (.1993.) . respectively. 



plots, we have ~ 300 stars in common. We can see a systematic 
offset between these two sets of results with a mean difference 
of -120 ± 5 (cr = 94) K. This offset is similar to that seen in 
ISousa et al ] (120081) . The comparison between the derived [Fe/H] 
seems to show closer agreement, although there is still a mean 
difference of -0.07 + 0.01 (cr = 0.13) dex that is significant. 

In the sam e figure, we also compa re our results with the 
photometry of ICasagrande et aP (l20I ih . In this case, to make 
a robust comparison of the two methods we only considered the 
stars in this work with more reliable photometry that are labeled 
with "irfm" pedigree (see Casagrande et al. (201 1) for more de- 
tails). In panel (c) of Figure |4] we can see the comparison be- 
tween the two temperatures for more than 200 stars. The results 



6. Solar twins in tlie sample 

The Sun is by far the most well-known star. But how similar 
are other stars to the Sun. To understand this question, a good 
starting point is to compare the Sun to similar stars. With this 
comparison, we can start to ask questions such as: How frequent 
are the stars in the solar neighborhood with the same mass, tem- 
perature, and metallicity? How different is their chemical abun- 
dance? Do they also host planets? Are these planets similar to 
the ones we find in our Solar System? 

T he quest to find solar tw ins has been ongoing for a long 
time (ICavreldeStrobelLll996 '). Here we present a few new solar- 
twin candidates within the volume-limited sample observed with 
HARPS. The list of star candidates to be solar twins are listed 
in Table [T] They were selected such that their "accuracy error" 
(given in Sect. 3) include the solar values for each parameter 
(5777K for effective temperature, 4.44 dex for logarithmic sur- 
face gravity, 0.00 dex for [Fe/H], and 1.0 for mass). From the 
table, one can see that the age estimates are also within the solar 
age. 
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Table 1. Solar twin candidates within the volume-limited sample observed with HARPS presented in this work. 



Star ID 


Tcfif 




log 






[Fe/H] 


N(Fe i,Fe ii) 


Age,, 


Mass,, 
[Mq] 




[K] 




[cms ] 


[km s"'] 








[Gyr] 


HD 19641 


5806 ± 


61 


4.39 


±0.10 


0.94 ± 0.01 


-0.01 


± 0.05 


255, 35 


4.6 ± 3.4 


0.992 ± 0.039 


HD 29263 


5780 ± 


65 


4.35 


±0.10 


0.94 ± 0.03 


0.03 


±0.05 


259, 35 


7.0 ±2.8 


0.995 ± 0.030 


HD 76440 


5764 ± 


63 


4.43 


±0.10 


0.89 ± 0.02 


-0.01 


± 0.05 


260, 35 


3.7 ± 3.2 


0.983 ± 0.038 


HD 96116 


5832 ± 


63 


4.52 


±0.10 


0.96 ± 0.03 


-0.01 


± 0.05 


260, 35 


2.7 ± 2.5 


1.005 ±0.036 


HD 134702 


5782 ± 


64 


4.50 


±0.11 


0.74 ± 0.04 


-0.04 


±0.05 


259, 36 


3.2 ± 3.0 


0.979 ± 0.037 


HD 145927 


5819 ± 


62 


4.41 


±0.10 


0.93 ± 0.02 


-0.03 


±0.05 


254, 35 


2.7 ± 2.5 


0.994 ± 0.036 


HD 216008 


5773 ± 


62 


4.38 


±0.10 


0.91 ± 0.02 


-0.04 


± 0.05 


257, 35 


6.1 ±3.9 


0.967 ± 0.037 



Notes. loggs,,ec the spectroscopic surface gravity; ^, is the microturbulance speed; N(Fei,Feii) is the number of lines used in the spectroscopic 
analysis; MasSp and Age,, are the mass and age determined directly from the Padova web interface using the Hipparcos parallax; The error 
presented here are the assumed "accuracy" errors, which were used to select solar twin candidates within the full sample. 
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Fig. 6. Comparison between the [Fe/H] derived from the 
CORALIE CCF and the [Fe/H] derived using our spectroscopic 
method. 



7. The largest volume-limited sample 

7.1. [Fe/H] based on the CCF calibration 

As akeady indicated in previous sections, the sample for which 
we present par ameters is an exte nsion of the CORALIE planet- 
search sample dUdrv et al[ 2000h. Thi s samp le was presented as 
a comparison sample in Santos et"an (l2004 l). where the metal- 
licities of all the stars were determined using a cross corre- 
lation function (CCF) calibration. We have already a spectro- 
scopic determination for a significant number of stars within the 
CORALIE sample for which we derived precise spectroscopic 
parameters using HARPS spectra analysed in previous work 



for other different planet search sa mples ISantos et al.^ |2004 
12005; "S ousa et al 112006'. 2008'. '2010V In Fig.|6l we compare our 
[Fe/H] spectroscopic measurements with the one derived using 
the CCF calibration for stars in the CORALIE sample. The re- 
sults shown in the figure for the 522 stars with a spectroscopi- 
cally determined metallicity are very consistent and have a mean 
difference of only 0.01 + 0.00 (cr = 0.07) dex. This small mean 
difference is negligible (when considering the typical [Fe/H] er- 
ror for each star) and its consistency allows us to use the CCF 
metallicities, for the stars without our spectroscopic derivation, 
to construct the largest uniform volume-limited sample by com- 
bining both the HARPS and CORALIE samples. This includes 
most FGK dwarfs in a volume of ~ 60 pc around the Sun. 



Although the consistency between our spectroscopic [Fe/H] 
and the one derived using the CCF calibration is very good, we 
use the more precise spectroscopic derivation for stars when one 
is available. This only concerns the CORALIE sample that we 
use here and is composed of a total of ~ 1250 stars for which 
522 stars have spectroscopically derived [Fe/H] in our previous 
works. When combining these stars with the HARPS sample 
presented in this work, we achieve a total number of stars of 
1830 (1 104 out of 1830 with homogeneous spectroscopic metal- 
licity measurements). This ensures that it is the largest sample of 
stars with precise and homogeneous [Fe/H] measurements ever 
studied in the framework of planet formation and evolution. 

7.2. Tlie frequency of planets with [Fe/H] 

In Fig. |2l we present metallicity distributions for the initial 
CORALIE sample (top panels), the HARPS samples discussed 
in this work (middle panels), and the large and combined sam- 
ple (lower panels). The middle and right plots correspond to the 
population of planet host stars in each sample. There are already 
a few planet hosts presented in the HARPS volume-limited sam- 
ple th at have been discov ered in the l ast years jT innev et af , 
'2003; Hebrard et al., 2010; Santos et al.', 2010; Naef et al., 201C; 
Moutou e t al.., 2009; L o Curto et al., 2010; Santos et al., ,20 11 ; 
Mordasi ni et akiboi U lMoutou e t 3111201 ll) . For all of the planet 
hosts, we were careful to include only the stars that have con- 
firmed planets. In this way, we obtained a total of 107 planet 
hosts, 80 in the CORALIE sample and 27 in the HARPS sam- 
ple. The information about the presence of planets was updated 
using the exoplanet encyclopedic^ and confirmed by the refer- 
ences therein. 

Looking at the top panels, for the CORALIE sample, the 
left panel shows the typical metallicity distribution for the so- 
lar neighborhood. This distribution is very similar to the ones 
presented in studies dealing specifically with the metallicity dis- 
tribution in the solar neighborhood which are based on unbiased, 
volu me-limited samples (e.g.Havwo od. 2001; Fuhrmann, 20041 
20081; ICasagrande et al.L [201 lb . The CORALIE distribution has 
a similar shape indicating that this sample is also rather unbiased 
and representative of the stellar population in the solar neighbor- 
hood. 

The distribution of the planet hosts is, as seen before, very 
different from that of the full sample, being more metal rich. 
The probability of belonging to the same distribution according 
to a Kolmogorov-Smirnov test (hereafter PKS) is nearly zero 



http ://exoplanet . eu/ 



6 



Sousa, S. G. et al.: Spectroscopic characterization of a volume limited sample. 





Fig. 7. Metallicity distributions for the CORALIE sample (top panels), the HARPS sample presented in this work (middle panels), 
and the union of the two samples (bottom panels). The left panels shows the full samples distribution of [Fe/H], while the center 
plots for each sample shows the planet host [Fe/H] distribution. The right plots shows the frequency of planets for each [Fe/H] bin. 



(PKS: 2.3e"^). The right plot shows the frequency of planet host 
for each bin of metallicity. This plots shows an evident higher 
probability of finding giant planets around metal-rich stars. 

The middle plots of Fig. |7] displays the same sequence of 
plots for the HARPS sample. This sample shows the same kind 
of the metallicity distribution as seen for the CORALIE sample. 
We find that both full samples have a very significant PKS of 
~ 83%. On the other hand, although the HARPS planet-host 
metallicity distribution is very different from that of the full 
HARPS sample (with a PKS of ~ 2%) it has a somewhat dif- 
ferent shape from the distributions in the top plot for the the 
CORALIE planet hosts. Their PKS is ~ 14%, a value that indi- 
cates that there are some differences between the two samples. 
However, when we use these data to plot the planet frequency 
for each metallicity bin within the sample (right panels of Fig. 
|7]) we can see that the plots are again quite similar, i.e. in both 
plots an increase in the planet frequency can clearly be seen for 
the metal-rich stars. 

Finally, in the bottom panels we see the same sequence of 
plots for both samples combined. The difference in the mean 
metallicities between the planet hosts and the full sample is 0. 15 
dex. The same trends observed before for each individual sample 
are also observed here, showing again the general correlation 
with metallicity that implies the giant planets are more likely to 
be discovered around metal-rich stars. 

The error bars expressed in the relative frequency plots were 
derived assuming a Gaussian distribution when the total number 
of stars in a bin was higher or equal to 100 or a Poisson distribu- 
tion when lower than 100 stars. 
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Fig. 8. Planetary mass distribution in log scale, illustrating the 
two evident peaks at two distinct planet mass regimes, the nep- 
tunian mass and the jupiter mass planets. 



7.3. Planet frequency, planetary mass, and [Fe/H] 

We analyze the planet frequency around stars considering the 
mass of the present planets. In Figure [8] we show the distribu- 
tion of the mass of the planets present around the stars in the 
large HARPS + CORALIE sample. The presence of two visi- 
ble peaks in different mass regimes is clearly visible, namely the 
Jovian planets and the Nept une-mass planets . This was previ- 
ously found in the review of iUdrv et al.l ( l2007h . In the same fig- 
ure, we show the intervals of mass of planets that we analyze. 
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Table 2. Number and frequency of stars and planet hosts for [Fe/H] bin considering the mass of the most massive planets present. 



[Fe/H] bin 


nstars 


np 


fnp (full) 


np 


fnp (A) 


np 


fnp (B) 


np 


fnp (C) 


np 


fnp (D) 


[-0.6, -0.5[ 


53 


2 


( 3.77 %) 


2 


( 3.77 %) 





( 0.00 %) 


2 


( 3.77 %) 





( 0.00 %) 


[-0.5,-0.4[ 


81 


1 


( 1.23 %) 


1 


( 1.23 %) 





( 0.00 %) 





( 0.00 %) 


1 


( 1.23 %) 


[-0.4, -0.3[ 


118 


2 


( 1.69%) 





( 0.00 %) 


2 


( 1.69 %) 





( 0.00 %) 





( 0.00 %) 


[-0.3,-0.2[ 


222 


6 


( 2.70 %) 


4 


( 1.80%) 


2 


( 0.90 %) 


2 


( 0.90 %) 


2 


( 0.90 %) 


[-0.2,-0.1[ 


296 


12 


( 4.05 %) 


11 


( 3.72 %) 


1 


( 0.34 %) 


5 


( 1.69 %) 


6 


( 2.03 %) 


[-0.1,+0.0[ 


336 


14 


(4.17 %) 


11 


( 3.27 %) 


3 


( 0.89 %) 


4 


( 1.19 %) 


7 


( 2.08 %) 


[+0.0, +0.1 [ 


291 


21 


( 7.22 %) 


21 


( 7.22 %) 





( 0.00 %) 


10 


( 3.44 %) 


11 


( 3.78 %) 


[+0.1,+0.2[ 


191 


17 


( 8.90 %) 


16 


( 8.38 %) 


1 


( 0.52 %) 


4 


( 2.09 %) 


12 


( 6.28 %) 


[+0.2, +0.3[ 


97 


21 


(21.65 %) 


20 


(20.62 %) 


1 


( 1.03 %) 


6 


(6.19 %) 


14 


(14.43 %) 


[+0.3,+0.4[ 


43 


11 


(25.58 %) 


11 


(25.58 %) 





( 0.00 %) 


4 


( 9.30 %) 


7 


(16.28 %) 


Total planets: 




107 




97 




10 




37 




60 




< FelH >: 
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Fig. 9. Metallicity distribution and respective frequency for the planet host with the four different planet-mass regimes described in 
Figure m In the left and bottom plot, we fit a power law to the data - see more details in the text. 



We consider four intervals of masses. The Jupiter-mass plan- 
ets with masses between 0.1 and 25 Jupiter masses (mass in- 
terval A). The Neptune-mass planets with masses between 0.01 
and 0.1 Jupiter masses (mass interval B). The Jovian planets are 
also sub-divided into the lighter jovians with masses between 0. 1 
and 1 Jupiter masses (mass interval C) and heavy jovians with 
masses between 1 and 25 Jupiter masses (mass interval D). 

In Figure |9] we present the metallicity distribution and the 
planet frequency for each of the mass intervals described. The 
frequencies of planets in each interval are also presented in Table 
|2] For this analysis and for the stars hosting planetary systems, 
we only consider the mass of the most massive planet present in 
the system. 

For interval A, which contains stars with giant planets we can 
again see the higher frequencies for the more metal.rich stars. 
In this case, we fit a power law to the histogram of the planet 
frequency. The fitted function is giving by: 

^{planet) 

which is not v er y diff erent from the one presented in 
IValenti & Fischeil (l2005h with the coeflicients (0.03, 2). 
Although in our case, the fitted power law does not appear to 
reproduce the very rapid increase in frequency around the 0.2 
[Fe/H] bin. This rapid increase is clearly evident, and may repre- 



sent a discontinuity in the planet formation efliciency as a func- 
tion of [Fe/H]. In the same figure, we can see the metallicity dis- 
tribution and planet frequency for the other intervals (B, C, and 
D). It is interesting to see a gradual increase in the mean metal- 
licity of each interval as we increase the mass of the planets. This 
may indicate that there is some superposition of the two groups 
of planets in the mass bins. We can also infer from this obser- 
vational result that the metallicity distribution of the Neptunian 
hosts is rather flat compared to that of the stars hosting jovians. 
This evidence illustrates once again the possible difference be- 
tween the formation of Nep tune-mass planets a nd J upiter-mass 
planets as already noted in ISousa et al] (l2008h .and lUdrv et aTl 
(2007), which is su pported by theoretical models based on the 
core accretion idea (IMordasini et al l l2009h . 

These results should be intrepreted with caution because ob- 
servational bias may be present at some level. On the one hand, 
the most successful methods such as the use of radial velocities 
or transits are biased in their detections because it is easier to 
find more massive planets and those closer to the host star. On 
the other hand, several planet-search programs are focused on 
specific samples, constrained to contain only high metal-content 
stars, or using observational strategies designed to identify short 
period planet. All these together can introduce observational bi- 
ases into the planet host subsamples. 
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8. Summary 

We have presented precise and homogeneous measurements of 
spectroscopic stellar parameters for a volume-hmited sample of 
stars observed within the HARPS planet-search program. The 
parameters were derived with a method similar to one we pre- 
viously used to derive a set of homogeneous and precise param- 
eters for all the stars. We estimated the mass of each star using 
the Padova web interface. 

We have verified the spectroscopically derived temperature 
using two other independent calibrations, one using the B-V 
color and measurement of [Fe/H] and the other using calibrated 
line-ratios. Both calibrations infer very consistent temperatures. 
We have compared our derived parameters with others found in 
the hterature. The comparisons performed here show compatible 
results, with some small offsets that were briefly discussed. 

We have combined the HARPS sample presented here with 
the original CORALIE sample. The latter sample is composed 
of stars that have either a spectroscopically determined [Fe/H] 
derived in our previous works, or a CCF-calibrated [Fe/H] that 
is proven to be compatible with our spectroscopic |Fe/H|. Both 
samples together constitute the largest volume-limited sample 
with homogeneous and precise [Fe/H]. 

Using the largest volume-limited metallicities, we have per- 
formed a simple analysis of the planet host distributions, focus- 
ing also on the masses of the planets present in each system. 
In this analysis we have confirmed again the higher frequency 
of planets for metal-richer stars, but we have also found a more 
rapid increase in frequency around 0.2 dex. 

After dividing the system accordingly to the mass of the most 
massive planet orbiting the planet host, we have found a clear 
trend with metalhcity. Stars hosting less massive planets have 
on average lower metallicity than stars with higher planet mass. 
This is clearly seen when dividing the system into three bins 
of masses: Neptune regime, fighter Jovians regime, and heavy 
Jovian regime. This result also should be interpreted with care 
because of possible observational bias present in the subsample 
containing the planet host. 
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